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Abstract 

The methodology of stochastic description for dissipation, a generic scheme to decouple the 
interaction between two subsystems, is applied to the study of dissipative dynamics in quantum 
optics. It is shown that the influence of the coupled thermal or vacuum field on the quantum mode 
can be exactly represented by the induced stochastic fields. The quantum mode thereby satisfies 
a stochastic differential equation and dissipation effect due to the coupling with the environment 
is obtained through statistical averaging. Within the framework of stochastic description, it is 
demonstrated how to derive the master equation for a single optical mode interacting with the 
bosonic bath. A numerical algorithm for solving the master equation in which the coefficients are 
determined by a set of integral equations is discussed and a comparison with the known results is 
displayed. The derivation of the master equation for the spontaneous decay of two-state atoms in 
the vacuum is also presented. 
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I. INTRODUCTION 



The best known dissipative dynamics is perhaps the Brownian motion, which has been 
greatly explored in theory and well understood since the pioneering work of Einstein 
Because the traditional Brownian particle is heavy and its surrounding environment or bath 
is at sufficiently high temperature, the motion of the particle can be accurately described 
by classical mechanics [2|. When the system of interest is very light, or the temperature 
of the environment is low, however, classical mechanics is no longer applicable and one has 
to invoke quantum mechanics. Actually, all physical systems intrinsically follow quantum 
mechanics and the traditional theory of Brownian motion should be a classical approximation 
of the exact quantum version. There have been many efforts made to establish a quantum 
formulation of dissipative dynamics and the most successful one is based on the system 
plus bath model [3]. Compared to classical counterpart, the quantum evolution exhibits a 
unique feature, that is, coherence. Quantum coherence is a consequence of the principle of 
linear superposition and plays indispensable role in the operation of quantum information 
devices |4|. It is also a fundamental issue related to quantum measurement 5KZ|. The main 
purpose of studies on quantum dissipation is thus to reveal how the environment affects the 
time evolution of the quantum system, in particular, the decoherence effect 

The classical Brownian motion is generally described by a generalized Langevin equation 
in which the random force is induced by the thermal bath 2|, |3J. One naturally wonders 
whether the impact of the bath can be defined as a classical random field. Kubo was the 
first to propose the stochastic Liouville equation for quantum dissipative systems, but his 
formulation is phenomeno logical js]. Given the system plus bath model, the key issue is then 
to acquire the equation of motion of the system, in which the dissipation effect due to the 
bath is exactly taken into account and the explicit dynamics of the bath does not show up. 
In other words, one aims at finding the effective motion of the system in its own physical 
space instead of inspecting the every detail of the dynamics for the whole system. To this 



end, several theoretical frameworks including the projection operator technique 



influence functional method 



12 



10 



, the 

14|, and quantum Langevin equation approach 15Hl7j] were 



suggested and applied to a wide range of physical systems. Of course, all of these methods 
have their own pros and cons in practice. We have recently put forward a stochastic theory 
for dissipative systems, in which the interaction with the bath is rigorously mapped into 



2 



stochastic fields during the evolution of the system 18| . In this formulation the system is 



181 ] . For comparison, in Kubo's 



subjected to complex stochastic fields induced by the bath 
stochastic Liouville equation, there is only a real stochastic field due to the bath. For spe- 
cific dissipation, we have proposed the hierarchical equation of motion approach in terms 
of the stochastic formulation 



19l . I20IL which has become an efficient, deterministic numer- 
ical technique of many applications [21| . Combining the stochastic and the deterministic 
methods, we were able to solve the dynamics of the two-state system strongly coupled to 
a bosonic bath [22]. Besides the trophy in numerical simulations, it has been shown that 
the stochastic formulation offers a convenient, systematic procedure for theoretical analysis, 
say the derivation of the master equation for linear systems 23j. When the existence of 
master equation is warranted, its derivation and solution should be the central task in quan- 
tum dissipative dynamics. This paper, as a continuation of the previous one [23j, addresses 
the dissipative dynamics in quantum optics. We will apply the stochastic formulation to 
derive the well-known master equations for two models. One is a single mode perturbed 



by a therma. 



vacuum 



18 



bath [24{ and the other is the spontaneous decay: the two-state atom in the 



25 



30]. 



As discussed in literature, the fluctuations of the thermal field are a major source of 
dissipation in quantum optics, which damage coherence of the system [31]. Again, because 
quantum optical devices operate when decoherence is negligible, to realize optimal function- 
ing, it is sometimes necessary to design a scheme to control dissipation 26, 31 -[33]. This 
entails a clear revelation of underlying dissipative dynamics. It is no doubt that exactly 
solvable models may provide fundamental understanding in this respect and are always de- 
sired. In the previous paper we have shown how to employ the stochastic description of 
quantum dissipation to work out the master equation [23] . The harmonic oscillator coupled 
to the Caldeira-Leggett heat bath is used as an example. In this paper the thermal field as 
well as the vacuum is considered to be the heat bath. There are two kinds of interaction 
between the system and the bath, one corresponding to the absorption and the other the 



emission, of a photon energy 24|, l34j-l36| . Although these models can formally be regarded 
as a result of rotating-wave approximation (RWA) imposed on the Caldeira-Leggett model, 
we will avoid the discussion on the validity of the approximation 

The first model we shall consider is essentially the d issip ative harmonic oscillator within 
RWA. Its master equation was derived by Xiong et al. 24(, resorting to the influence func- 
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tional approach developed by Feynman and Vernon [3, Q]- As the authors observed, when 
the coupling between the system and the bath becomes strong, the dissipative dynamics will 
change dramatically because of the non-Markovian memory effect 2^]. In a recent paper jssj, 
Tan and Zhang used the same method to discuss the consequence of initial system-bath cor- 
relation. The second model describes the spontaneous decay of two-state atoms in vacuum, 
which is exactly solvable. It has been frequently used as benchmark calculations in quantum 
optics. Indeed, this model is so well-known that diversified methods have been used to derive 
its master equation 



18 



25 



30]. 



The paper is organized as follows. In Sec. II we recapitulate the stochastic formulation 
for quantum dissipation. In Sec. Ill we apply the scheme to a single-mode cavity system 
coupled to a thermal field and derive its exact master equation. In Sec. IV the obtained 
master equation is shown to be equivalent to the result in Ref. [24j and some remarks on 
the numerical implementation are given. In Sec. V the master equation of the same system 
subjected to a driving external field is derived. In Sec. VI we show how to derive the master 
equation of a two-state atom coupled to the vacuum field. We present our conclusions in 
Sec. VII. 



II. THEORY 

To study the dissipative dynamics of a quantum mode in an optical cavity, we start 
with an arbitrary cavity system coupled to a thermal field consisting of infinite number of 
harmonic oscillators. The Hamiltonian of the entire system assumes 

h = h s + J2 M^- + E h { c M + c .M) , (i) 

3 3 

where the first term on the right-hand side is the Hamiltonian of the cavity mode, the second 
term is the Hamiltonian of the thermal bath, and the last two terms define the interaction 
between the system and the bath. Here f\ and are operators for the system and they 
are a hermitian pair, f 1 = Note that the two interaction terms can be interpreted as 
emitting and absorbing a quantum phonon or photon by the bath. The model will be the 
Caldeira-Leggett type when the interaction is of the form £\ hcj (fx + fi j (b^ + bj\ . As 
shown in the previous papers [3], the dissipative dynamics can be described by a stochastic 
formulation in which the system evolves in the stochastic fields induced by the bath and the 

4 



statistical average of the random density matrix is nothing but the reduced density matrix. 
For the model we consider, the random density matrix satisfies 



ihdp s {t) = 



H s + ^2gk(t)fk,p s 



k=i 



dt + 



£[A, 



Ps 



k=l 



dW lk + i^-^2{^P s } dW ^ ( 2 ) 

k=l 



where the bath-induced stochastic fields are given by 

Tr b {h Cj b]p b (t)} 



9i 



(3) 



and 



92 



w = E 



Tr b {hcjbjp b (t)} 
Tr b {p b (t)} 



(4) 



Here, introduced are the complex Wiener processes W lk {t) = f*dt'[v lk (t') + iv^{t')\ and 
W2k(t) = Jq dt' \y 2k {t') + iv3k(t')]i where v nk (t) (n — 1 — 4) are independent Gaussian white 
noises with zero mean and delta function correlation. This is the main result of the stochastic 
formulation and will be the working formula. To use it, of course, we need to first calcu- 
late g~i(t) and g 2 (t). In this formulation, gi(t) and g 2 (t) can be determined by the evolution 
of the bath, 

hy/hsr^ r,+ 1 j h\ffi 



ihdp b = ^2 [Hujtfjbj,^ dt + — ^— ^ c j [ b ]>Pb\ dW 2 i + — ^— E Cj [bj,Pb] dW 22 



The formal solution of p b (t) can be written as 

p b (t)=u 1 (t,0)p b (0)u 2 (p,t), 
where ui i2 (t,0) are the forward and backward propagators dictated by 

^(*) = E + E c j 6 ]^i(*) + E c j 6 ^i2(o 



(5) 



(6) 



and 



h 2 (t) = 2^ ^jbjbj + 2^ CjbjV2i{t) + ~y~ 1_j c i b i 

j 



iV22(t). 
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with 



7/H (t) = U Z l{t) + iv 3 x(t) + Wu(t) + 1/41 (t), 

^L2(*) = ^22 (*) + iVaait) + ^12 (*) + ^42 (*), 

7/21 (t) = f2l(*) + (0 - - ^41 (*), 

7/22 (t) = I/ 2 2 (t) + ^32 (0 ~ ^12 (t) ~ ^42 (*) , 

being complex white noises. Because the bath modes are independent, the propagator of 
the bath is a product of the individual ones, namely, Ui(t, 0) = Ylj M j,i(^0) and «2(0,t) = 

rii«i,2(o,t). 



As illustrated in the previous paper and other references [23j, |39| , the propagator for each 
bath mode can feasibly be obtained upon using the interaction representation. As a result, 
the forward propagator Uij(t, 0) reads 



u jA (t,0) = C i|W (t)e cr * ll ^e * ia « 6 J« il o(t > 0) ) 



(7) 



where Uj 7 o(t, 0) is the propagator o 



hujjb^bj, which is well known [l2 
Cj,io(0 = ex P 



13 



the undriven harmonic oscillator described by hj t o 



40 



43|, and 



* Jo 

* Jo 

Similarly, the backward propagator Ujp(0, t) is 



(8) 



where 



~Z C ?/ dtl j dt ^2l(tl)V22(t2)e^ 



Cj,2o{t) = exp 

* Jo 

C h22 {t)=i^c 3 fdhrmWe-^-V. 
z Jo 
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Assume that the thermal field is initially in thermal equilibrium, 

Mo) = V"4lI e " %0 . ( 9 ) 

3 

where Z b = Ti h is the partition function. Inserting together with Eqs. (7) and (8) 

into Eq. (6), carrying some operator algebra and rearranging, we obtain 

Pb(t) = II F i® ex P { [°iMt) + CiMt)e-^] b]) exp { [C jA1 (t) + C j>21 (t)e^] 6, } exp {-(5H b ) , 

3 

(10) 

where 

F 3 it) = ^C hW {t)C h2Q {t) exp {C,- u (f) [C hl2 (t) + C m {t)e~^\ } . 
Then the bath-induced fields defined by Eqs. (3) and (4) can be worked out, 

9i(t) = ^ f dt' { ai (t - o MO + 11*2(0] -a 2 (t- [^(0 + MOD (ii) 

and 

= jT - 0M0 + ^3i(01 + - OKi(0 + MO]}, (12) 

where ai(t) and a^) are response functions given by 

roc 

ai (t) = / dcuJ(co)e iuJt (13) 
Jo 

and 

a 2 (t) = jf duJ(tu) coth e^, (14) 

J(oj) being the spectral density function 

J( W )^cj%- W ). (15) 

We like to stress that J(w) completely captures the influence of the thermal field on the 
system. With the known gi, 2 (t), Eq. (2) is a closed stochastic differential equation. That is, 
one can generate the required stochastic fields through white noises and solve Eq. (2) for a 
given initial condition p s (0). The reduced density matrix p s (t) is of course the mathematical 
expectation of random density matrix p s (t), namely, p s (t) = M {p s (t)}. We can also try to 



find the equation of motion for p s (t) from Eq. (2). To this end, we take stochastic averaging 
to obtain 



dt 



H S} p 
ihy/h 



+ 



ihyh 



fx, / dt' a x {t - t')6 s>11 {t,t') - a 2 {t - t')6 S)12 {t,t') 



f 2 , I dt' al(t-t')0 Sj21 (t,t') + a*(t-t')0 Sj22 (t,t') 
o 



(16) 

(17) 
(18) 
(19) 
(20) 

In the above derivation the nonanticipating property of Ps(£L namely, 
M {ps(t)dW mt k(t)} = (m = 1,2), and the Furutsu-Novikov theorem [43], that is, 
M {y(t')F [u]} = M {8F[u]/8u(t')} for a white noise v(t) and its arbitrary functional F[u], 
are used. As shown in Ref. [18j], the formal solution of Liouville-like Eq. (2) can be written 

as 

Ps (t) = fA(t,o)p s (o)f/ 2 (o,t), 

where Ui(t,0) is the forward propagator associated with the Hamiltonian 



where the dissipative operators are 






a,n(M') 


=M < 


f 5 Ps {t) 5 Ps (t) 
[5is 22 (t') "5v 32 (t>) 


O s , 12 (t,t') 


=M< 


r, s Ps (t) s Ps (t) 

[ 5u 12 (f) 8v 42 (t>) 


4,21 (M') 


=M< 


f 5 Ps (t) 5 Ps (t) 
|> 21 (tf) "Su 31 (f) 


6 s , 22 {t,t') 


=M< 


\.5p s (t) 5 Ps (t) 
[Vi(t') Su 41 {f) 



(21) 



H 1 (t)=H a + £ 11 (t)h + {v i (t)f 2 
while U 2 (0,t) is the backward propagator associated with the Hamiltonian 

H 2 (t)=H a + Z2i{t)f 1 + Z 22 {t)f 2 

with 



(22) 



(23) 





= 9x{t) 






= Ht) 




6i(0 


= 9x{t) 






= 92® 
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Following Ref 18j , we calculate the functional derivatives to obtain formal solutions of the 



dissipative operators, 

a,n(M') = 4=M [u 1 (t,t')f 2 U 1 (t', t)p 8 (t) + p a (t)U 2 (t, t')f 2 U 2 (t', *)}, (24) 

O s , 12 (M') = {u^t^hU^t'rfpsit) - p 8 {t)U 2 {t, t')f 2 U 2 (t',t)}, (25) 

6 S!21 (t,to = 4^ M (^(* ) *OAt / i(^*)p s (*) + p s (*)t / 2(t ) t / )A^(t',*)), (26) 



6 S , 22 (M') = -^M {[/!(*, -p s (t)U 2 (t,t')hU 2 {t',t)} . (27) 

When these operators can be expressed in terms of the reduced density operator p s (t) and 
other operators of the system, then Eq. (16) becomes a master equation. In the following 
section, we show that a dissipative single-mode optical cavity is indeed the case. 

III. MASTER EQUATION FOR DISSIPATIVE SINGLE-MODE OPTICAL CAV- 
ITY 

Let us consider a single-mode cavity, H s = hw a^a, with the coupling operators fx = 
a and f 2 = a) . Therefore, the dynamics of the random cavity are determined by the forward 
and backward Hamiltonians 

Hx{t) = hcoo^a + £ n (t)a + Z 12 (t)<J (28) 

and 

H 2 {t) = hw Q a ] a + £ji(t)a + fa(t)a*. (29) 
These Hamiltonians are exactly solvable and one readily finds 

Ux(t,t')aUx(t',t) =ae ib *W + % - f dt^^e^-^, 



U 2 {t,t')aU 2 {t',t) =ae**>W + % - \ dtx^^e^ ^, 



U^t'^Uxit' ,t) =a^e- lwo{t - t ' ) - % - j dti6i(*i)e" <Wo( * 1 "*' ) , 



U 2 (t,t')a^U 2 (t',t)=a^e~^ t - t ^ - % - \ dtx^i^e'^-^ . 



Inserting into Eqs. (24) — (27) and taking statistical averaging, we obtain 

6 s ,n(M') =^= W,p s (t)} e -^(*-*') + jf*^ jT di 2 e- <w °^- t ')a 1 (i 1 - * 2 )O a>11 (*,* 2 ) 

- /V /W-™ (tl -*W*i-*2)d ai i2(M 2 ), (30) 

O a ,i 2 (*,O=^[a t ,P.(0]e- iW0(t -* ,) -^ rfti^" ^e-^^-^aKta-fOa.ia^fa), (31) 

6 s , 21 (t, t') =4= {a, &(*)} e*" ^ + r d*i f 1 dhe^-^aKh - t 2 )6 8 , 21 (t, t 2 ) 
yn Jf Jo 

+ f dh f dhe^-^^ih-tJOs^t,^, (32) 
Jf Jo 

O a , 22 (*,O=4=k/5.(0]e <W0(t - t ' ) - / <**i / d* 2 e <W0(tl - t ' ) ai(t2-ti)O 8) 22(t,t2). (33) 
V« ./i' Jti 

In the above derivation, the following functional derivatives with respective white noises are 

used, 

M {&) + = - in t dhai{t ' - e)6 '^ 

M + } = " m t M - t>)d - M 
M { i iS + ll)} =ifi /* Ql(tl - 0<5 - (Ml) ' 

which can be found directly through the solution of p s (t). Note that O sll (t,t') = 
Sj 2i(t,t') and O sl2 (t,t') = —O s>22 (t,t'). By iteration, one can show that the dissipative 
operators should assume the following forms, 

d s ,ii(M') =zn(M') {a\p s (t)}+x 12 (t,t') [a\p s (t)] , (34) 

O s , 12 (t,t') =x 21 (t,t') [a\p s (t)] , (35) 

O s , 21 (t, t') =x* u (t, t') {a, p s (t)} - x* 12 (t, f) [a, p s (t)\ , (36) 

O s ^(t,t')=x* 21 (t,t')[a,p s (t)]. (37) 

Because the operators {a^,p s (t)}, [a*, p s (t)] , {a,p s (t)}, and [a, p s (t)] are arbitrary, it 
is straightforward to find out the equation of motion for x m k(t,t') (m,k = 1,2) with 
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Eqs. (30) — (33). The results are 



%(M') =^ e -^t-t>) + f dh / 1 dt 2 e--o(*i-*') Q;i ( fl _ f2 ) Xll (t )f2 ) ) (3 8 ) 

Vh Jt> Jo 



(39) 
(40) 



x 12 (t,t')= I dt x \ d^e-'^-^a^ -t 2 )x 12 (t,t 2 ) 
t' Jo 

dh / dt 2 e- iuoi - tl - t ' ) a 2 {t 1 -t 2 )x 21 {t 1 t 2 ) 1 



x 21 (t,t')=4=e^ o( ^' ) - / dt x [ dt 2 e~ iuo( *-V a *(t 2 -t 1 )x 21 (t,t 2 ). 
Vh Jt> Jt! 



Whenever these coefficients are solved, the dissipative operators Eqs. (34) — (37) become 
available. Inserting into Eq. (16), we obtain the master equation, 

dp s (t) 



dt 



- iAi(t) \p}a, p s (t)] + A 2 (t) [2ap s (t)a' - a*ap s (t) - p s (t)a^a] 
+ A 3 (t) [a'p s (t)a + ap s (t)a t - aJap s (t) - p s (t)aa t ] , 



(41) 



where the coefficients Aj(t) (j = 1 — 3) are defined by 



A^t) =u + lm 



A 2 (t) =Re 



A 3 (t) =Re 



Vh [ dt'al(t-t')x* n (t,t') 
Jo 

Vh / dt'a\{t-t')x* n {t,t') , 
Jo . 

Vh i dt'a* 2 (t-t')x* 21 {t,t') -Re Vh dt'a\{t - t')y*(t,t' 
Jo _ . Jo 



(42) 
(43) 
(44) 



with 



y*(t,t')=x* n (t,t')+x* 12 (t,t'). (45) 

It is clear that Ai(t) is a frequency- renormalization coefficient, A 2 (t) and A 3 (t) are related 
to the conventional dissipation and fluctuation coefficients, respectively. In the following 
section, we will show the equivalence between our derived master equation and that by 
Xiong et al. in terms of path integral approach 24]. 



IV. COMPARISON WITH KNOWN RESULTS 

Resorting to the influence functional method developed by Feynman and Vernon, Xiong 
et al. elaborated the derivation of the master equation of optical cavity coupled to a heat 



bath 



241 ] . For the case of a dissipative single-mode, their result is of the same form as 
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Eq. (41) and the corresponding coefficients read 



£i(t) = Wo + Im 
B 2 (t) = Re 
B 3 (t) = Re 



Vh [ dl/al(t-t')xu(t,lf) 
Jo 

Vh / dt'al{t-tf)x n {t,t') 
Jo 

Vh I dt'a* 2 {t-t')x 21 {t,t') 
Jo 



Vh [ dt'a\{t -t')y(t,t') 
Jo 



(46) 
(47) 
(48) 



The functions x n (t,t'), x 21 (t,t'), and y(t,t') are defined by 

%(M') = -^=u{t')u-\t), 
x 21 (t,t') = -±= u *(t-t'), 
y(t,t') = -L«*(0 - A [ U (t')u-\t)v(t) - v(t')} , 
where u(t) and v(t) obey the following integro-differential equations, 



u(t) + ico u(T) + / rft'a*(r - t')u(t') = 



(49) 
(50) 
(51) 



(52) 



and 



v(r)+iou v(T)+ [ dt'al(r -t')v(t') = \ [ dt' [a* 2 {r - t') - «*(r - t')} u* (*') (53) 
Jo 2 Jo 

with the initial conditions -u(O) = 1, v (0) = 0, and -u(r) = -u(t — r). To prove the equiva- 
lence of the results obtained by two different methods we only need to prove that Aj(t) = 
Bj{t) (j = 1-3), respectively. As displayed in Eqs. (42) -(44) and Eqs. (46) -(48), all 
definite integrals in the functions Aj(t) and Bj(t) are taken over the same time range [0,t\. 
Therefore, a sufficient condition for Aj(t) = Bj(t) is that the corresponding integrands are 
identical. Moreover, because these integrands consist of the factors «i(t) and a 2 (t) that are 
dependent on the specificity of the dissipation and can be arbitrary, one can further simplify 
the problem as a proof of following relations, 



x* u (t,t') = x n (t,lf), 
x* 2l (t,t') = x 21 (t,t'), 
y*(t,t') = y(t,t'). 



(54) 
(55) 
(56) 
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A. Proof of A 1 (t) = Bx (t) , A 2 (t) = B 2 (t) 



As clarified above, if Eq. (54) holds, then A^t) = B^t), A 2 (t) = B 2 (t). Note that u(t) sat- 
isfies the linear integro-differential equation (52) and that x n (t,t') = u(t') / (\fhu(t)). When 
the first argument t is fixed, xn(t, t') can be seen as a function of the time variable t', which 
also obeys Eq. (52), namely, 

d /■*' 

— x n (t, + iw xu(t, f) + J dhaKt' - h)x(t, h) = 0. (57) 

Return to the integral equation of x n (t,t'), Eq. (38). Calculating the time derivative with 
respect to t' and taking the operation of complex conjugation on both sides of Eq. (38), 
one obtains for xl^tjt') the same equation as Eq. (57). Also, the initial condition for 
these equations are the same, namely, x^t, t')\ t / =t = x n (t,t')\ t i =t = l/y/h. Therefore, 
A^t) = B^t) and A 2 (t) = B 2 (t) are proved. 



B. Proof of A 3 (t) = B 3 (t) 

One only needs to demonstrate that Eqs. (55) and (56) hold. A straightforward algebra 
shows that x 21 (t,t') is time-translation invariant, i.e., x 21 (t,t') = x 2 i(t + A, t' + A), where A is 
a constant. It means that x 2 \ is a function of the time difference t — t', x 2 i(t,t') = x 2 i(t — t'). 
As a result, Eq. (40) can be simplified as 

x 21 (s) = -±=e- i0J0S - [ dh f dt 2 e- i0JQtl al(t 2 - tjx^s - t 2 ). (58) 

\h Jo Jti 

Taking the first-order derivation with respect to the argument s and the complex conjuga- 
tion, one obtains 

d f s 

-j^x* 21 (s) = iu x* 21 (s) - J dt-yaKti - s)x* 21 {t l ), (59) 

subjected to the initial condition X2i(s)| s =o = By definition Eq. (50), the func- 

tion x 2 i(t,t') is only dependent on the time difference s — t — t' . Taking the operation of 
complex conjugation on both sides of Eq. (52) leads to the equation which is the same as 
Eq. (59). Besides, x 2 \(s)\ s= o = l/y/h = X2i( s )|s=o- Therefore, x 21 (t,t') = x 21 (t,t') does 
hold. 

By definition Eq. (51) and with the help of Eqs. (52) and (53), we find the that y(t, t') sat- 
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isfies 



d_ 

at 



7 y(t,t') = -iu y(t,t') - [ dtxal(t' -t 1 )y(t,t 1 )+ [ dt x al{t' - ti)5 2 i(Mi). (60) 

Jo Jo 



The same equation can be obtained for y*(t,t') from Eqs. (38) and (39). Moreover, 
y(t, t')\ t i =t = y*(t,t')\ t i =t = l/\/h. Therefore, one proves y(t,t') = y*(t,t') and as a re- 
sult, A$(t) = Bs(t). We have therefore demonstrated that the master equation Eq. (41) 
resulting from stochastic description is identical with that derived with influence functional 



method 



24J. 



Some remarks on the calculation of the coefficients of the master equation are in order. 



As discussed above, our procedure provides a set of integral equations, while Xiong et al. [24 ] 
introduce an integro-differential equation or the equation of the related Green's function. 
It is straightforward to numerically independent both of the two schemes to determine the 
coefficients. Although we prove that these two frameworks give the identical results, their 
numerical performance might be different. Because the computational scaling for solving 
the integral equation is less favorable than solving the corresponding differential equation, 
the Green's method is preferred in practice. 



V. DRIVEN CAVITY DYNAMICS 



Let us consider the cavity dynamics in the presence of a time-dependent external field e(t). 
Now the Hamiltonian of the system reads H s (t) = huoa^a + e(t) (a + a^). The master 
equation can be derived along the same line discussed in Sec. III. Although the external field 
only directly acts on the cavity system, and does not change the bath-induced stochastic 
fields, it does interfere with the bath during the evolution of the system. This effect is 
reflected in the change of dissipative operators. Starting with Eqs. (24) — (27), we solve the 
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required propagators and take the stochastic averaging to obtain 

4,n(M') =^= W,p s (t)} e-^-^ + J\h £ dhe'^-^a^h - t 2 )6 8 , 11 (t,t 2 ) 

Jt> Jo nvnJt' 
6 8 , 12 (t,t') =^= [a\p 8 (t)] e --°(*-*') - jf d*! jf diae-^^-^aKia - fi)d M2 (M 2 ), 



'ti 

6 S , 21 (M') =-^= {a,p s (t)}e^-<') + jf dix jf 1 dfae^-^aK*! - t 2 )d a , 21 (t,t 2 ) 

+ f dt x [\t 2 e i "°^a 2 (t 2 -t 1 )d 8 , 22 (t,t 2 ) + ^ f dhe^-^eih)^), 
Jf Jo nynJv 

O s , 22 (t,t') =^= [a,p.(t)] e^-V - jf d tl jf d^e^^-^a^^ - ti)6 8)22 (t, t 2 ). 

We use the same reasoning as that in Sec. Ill to obtain 

d s , 11 (t,t')=x 11 (t,t'){a\p s (t)}+x 12 (t,t') [a\p s (t)} + x 13 (t,t')p 8 (t), 
O s , 12 {t,t') = x 21 (t,t') [a f ,p s (0] , 

a,2i(i, = x* u (t, t') {a, p s (t)} - x* 12 (t, f) [a, p s (t)} + x* 13 (t, t')p 8 (t), 
O s , 22 (t,t') =x* 21 (t,t) [a,p s (t)\, 

where all coefficients except Xi 3 (t, t') are the same as that of the undriven case 
[Eqs. (38) — (40)]. The additional new function is defined by 

x 13 (t,t') = -- ^= / dhe-^-t'hih) + / dh ^ dhe-^-^^ih-t^x^h), 

frV h Jv Jf Jo 

which is linearly dependent on the external driving field e(t). 
With these expressions the master equation now reads 

= [-iA^t^a + C(t)a + D(t)a\ p 8 (t)] + A 2 (t) [2ap s (t)a^ - a!ap 8 (t) - p 8 (t)a!a] 
+ A 3 (t) [a ] p s {t)a + ap s {t)a ] - a ] ap s {t) - p s (t)cuJ] , (61) 



where 



and 



C(t) = - % -e{t) + dt' ai (t - t')x 13 (t, f) 

D(t) = --e{t) - ^ J*dt?al(t - t')x* 13 (t,t'). 
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Here, the coefficients Aj(t) (J = 1 — 3) are the same as that of the undriven case, which 
satisfy Eqs. (42) — (44). It becomes clear that there are effects of the external field on the 
system, one is the direct interaction and the other results in the very interplay between the 
driving field and dissipation. As a consequence, the external field can be applied to control 
dissipation, or via versa, dissipation can be used to modulate the external field. 



VI. MASTER EQUATION FOR TWO-STATE ATOMS IN VACUUM 

The spontaneous decay of a two-state atom coupled to a vacuum is described by the 
Hamiltonian Eq. (1) with H s = —hwoa z /2, fx = a~, and / 2 = a + , where a z is the pauli 
matrix, and a + and a~ are the raising and lowering operators. They satisfy the commutation 
relations [a + ,cr _ ] = a z , [a + ,a z ] = —2o~ + , and [cr~,o" 2 ] = 2<r~. This damped two-state 
model might provide fundamental understanding of decoherence and other features of the 
dynamics of a qubit coupled to a heat bath. It is no wonder that its master equation has been 
derived and explored by several authors with diversified theoretical methods. For instance, 
Garraway developed a pseudomode technique to solve the dynamics [25]. Through the 
solution of the Schrodinger equation for the entire system, Breuer and coworkers worked out 
the reduced density matrix and thereby proposed a simple derivation of the corresponding 
master equation by a brute force calculation of the derivative with respect to time 26|, 27]. 
They also developed a stochastic wave function approach to simulate the dynamics 28 ]. 



Strunz et al. proposed a different stochastic Schrodinger function method to solve the 



dissipative dynamics of the model 29(. The reduced density matrix resulting from the 
Schrodinger equation was also exploited by Vacchini and coworkers who recently showed 
how to generate the exact master equations corresponding to the time-convolutionless form 



and to the Nakajima-Zwanzig non-Markovian form 30|. In the first paper on the stochastic 



description of quantum dissipative systems, one of the authors a 
derive the master equation from the stochastic equation of motion 



so demonstrated how to 



18| . His method is based 



on self-consistency of an ansatz related to a stochastic average and the derivation was not 



expounded in the paper 



18| 



It seems that all the derivations in the literature are not direct and straightforward 
within one theoretical framework. We will show the stochastic description does offer a 
good pass to the master equation from the equation for the random density matrix for 
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the system. Because the bath is the vacuum field, the temperature is zero. As a re- 
sult, coth [fku I \2k B T)} = I and the response functions derived by Eqs. (13) and (14) become 
identical, cti(t) = a 2 (t) = a(t). Therefore, the bath-induced stochastic fields determined by 
Eqs. (11) and (12) become 



9i(t) 



dt'a(t - t') [-^12 (0 + 1^(0 + 11/32(0 - ^42 (£')] 



and 



Ut) = P J dt'a*(t - [iv u (t') + v 21 (t>) + iv 31 (t') + !/« (0] • 



(62) 



(63) 



The formal solution of the random density matrix of the system is still given by Eq. (21) 
where the forward and backward propagators U\(t,0) and 1/2(0, t) are ruled by the corre- 
sponding Hamiltonians Eqs. (22) and (23) with /1 = a~ and f 2 = <j + . Apparently, there 
are six complex Gaussian fields, g~i(t), g 2 (t), ^n(t), r}* 2 (t), rj 21 (t), and rj 22 (t) involving in the 
dynamics. Note that all of the six Gaussian noises have zero means and null autocovari- 
ances. The average of a stochastic process generated, therefore, is fully determined by their 
non- vanishing cross-covariances. Given g~i(t) and g 2 (t) by Eqs. (62) and (63), however, one 
can readily check that the white noises r]* 2 (t) and r] 21 (t) are not correlated with other four 
and do not have any influence on the averaged dynamics. Therefore, r)l 2 (t) and r] 2l (t) can 
be safely omitted when calculating the reduced density matrix. 

To derive the master equation, we insert g~i(t) and g 2 (t) into Eq. (2) and take stochastic 
averaging to obtain 



dm 
m at 



H s ,p s (t) 



+ 



2 

t 



a 



a + , [ dt' 
Jo 



dt'a{t-t')O sA {t,t') 
a*{t-t')6 s , 2 {t,t') 



(64) 



where the dissipative operators are 



6 Syl {t,t') = M 
2 

~ 71 



S Ps {t) 6p 9 (t) , . 6p 9 (t) 5 Pa (t) 



Su 12 (t') 5v 22 {t>) 
M{p s (t)at(t,t')} 



+ 1 



Su 32 (t') Su 42 (t') 



(65) 



and 



O s , 2 (t,t') = M\t 



Sp s {t) S Ps (t) , . 5p s (t) ^ 5 Ps {t) 



Svu(t') Su 21 (t') 5v 31 (t>) Su 41 (t') 



(66) 
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with af 2 {t,t') = Ui :2 (t, t')a ± Ui t2 (t' , t). We like to stress that the derivation up to now 
is parallel to that illuminated in Sec. III. Now we need to find the explicit expressions 
of M {p s (t)o'£(t,t')} and M {erf (t, t')p s (t)} in terms of p s (t) and other known operators of 
the system. To this end, we consider their derivatives with respect to t' , 

{p s (t)at(t,t')} = - iu M {p s (t)a+(t,t')} 

+ [ dha{t' -h)M {p s {t)at{t,h)U2{t,t')(T z U2{t',t)} (67) 
Jo 

and 

^M{a^(t,t')p s (t)} =iu M{ai(t,t')p s (t)} 

+ / dhofit' -h)M {U^t^aMt'^it^Psit)} ■ (68) 
Jo 

By virtue of a z = 2<7 + a~ — I, the two equations can be converted to 

{p s (t)a+(t, t')}=- iu M { Ps (t)a+(t, t')} + 2 J dha(t' - h)M {^(t, h,t') } 

- / dhait' ~h)M {p s {t)at{tM)} (69) 
Jo 

and 

{<r-(t, t') Ps (t)} =iu M {<r-(t, t')p s (t)} + 2 J dha*(t' - h)M {x 2 (t, t^t')} 

- [ dt ia *(t' -h)M {crf(Mi)p s (0} , (70) 
Jo 

where 

Xi(Mi,0 = p s {t)at{tM)^t') 

and 

x 2 (t,t u t') = ^i(t,o^r(*>*0^(*), 

with a lt2 (t,f) = U lt2 (t,t')a + a-U lt 2(t',t). 

To find closed equations for M {p s (t)cr^(t, £')} and M {erf (t, t')p s (t)}, therefore, we 
should evaluate M j.Xi(t, ti, t')| and M |x 2 (t, ti, t')|. When the first argument t is fixed, 
M |Xi(t, ti, £') j and M |x 2 (t, £ l5 t')| can be taken as the functions of t\ and t'. For brevity, 
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the argument t for functions X x and X 2 will not be written. On taking their derivatives 
with respect to t 1 , we obtain 



_d 
dt 



-M{Xx (*!,*')} =-iw M{x 1 (t 1 ,O} + y dt 2 a(*i-*2)M{2X! (t 2 ,*i)^ 2 (*,0 -*i(*2,0 

(71) 



and 



9t 



(72) 



We like to point out that the solutions for X\{ti,t') and X 2 (ti,t') can be many as long as 
their stochastic averages satisfy Eqs. (71) and (72). Because any solutions are sufficient for 
our purpose, we only consider the simple ones determined by 

l-X^t') = -iu X 1 (t 1 ,t') + [ dt 2 a(h-t 2 ) hx^MfaW) -M^t')} (73) 
oti Jo L J 

and 

^(ii/) =iu X 2 (t 1 ,t) + j dt^ih-h) [2ai(t, t')X 2 {t 2 M) - X 2 (t 2 ,t')j , (74) 

with the initial conditions Xi(ti,t')\ tl =t' = and X 2 (ti,t')\ tl =t> = 0. 

As a result, we obtain X^t^t') = and X 2 (ti,t') = 0. Then Eqs. (69) and (70) become 

d_ 

dt 

and 



-M {p s (t)at(t,t')} = -too M{p s (t)a+(t,t')} - / dt l0i (t' - h)M {p a (t)<T+(t, h)} (75) 



— M{o-i(t,t') Ps (t)} =tu M {o-^{t,t')p s (t)} -J dha*(t' -h)M {a^(t,h)p s (t)} . (76) 
They are integrated over time t', namely, 

M{p s (t)at(t,t')} =e- l ^ t '-^~p s {t)a + + f dt, P dhe'^' ~^ a(h - t 2 )M {p s (t)at(t,t 2 )} 

Jt 1 Jo 



(77) 



and 



M{a^(t,t')p s (t)}=e iuo ^a-p a (t)+ [ dt, I * dt^' ~^ a* {h - t 2 )M {a^(t,t 2 )p s (t)} . 

Jt' Jo 

(78) 
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Note that M {p s (t)a^(t,t')} j = M {crf(£, t')p s (t)}. By iteration, one can find 
that M {p s (t)a£(t,t')} and M {a{ (t, t')p s (t)} posses the following forms, 



M {p s (t)a+(t, t')} = x(t, t')p s (t)a 



(79) 



and 



M{ar(M> s (£)} =x*{t,t')a-p 8 {t). (80) 

Because the operators p s (t)a + and a~ p s (t) are arbitrary, the coefficient x(t, t') is determined 
by Eq. (77), which obeys 



t ft 



x (t,t') = e-^t'-t) + / dh / dhe-^'-^aih-hjxit,^). (81) 

Jt' Jo 

With the explicit expressions of M {p s {t)a+(t,t')} and M {a^{t, f)p 8 {t)}, Eq. (64) imme- 
diately becomes the resulting master equation. For the spontaneous decay of a two-state 
atom it reads 



dp s {t) 
dt 



i 

h 



H s ,p s (t)] [a + a-,p s (t)] 



+ R(t) 



a Ps{t)o + - ifj+a p s {t) - ]^p s (t)e + o~ 



(82) 



where S(t) and R(t) are the time-dependent coefficients for the descriptions of a frequency 
shift and a decay rate, respectively. Their expressions are 

ft 



R(t) =2Re 



dt'a*(t-t')x*{t,t') 



and 



S(t) =2Im 



dt'a*(t-t')x*(t,t') 



VII. CONCLUSION 



The main goal of investigating dissipative systems is to solve their properties, in par- 
ticular, to reveal the dissipative effect on their dynamics or Brownian motion. From the 
system plus environment model, we have shown 18j that the coupling to the environment 
can be rigorously mapped into stochastic fields and thereby provided a microscopic descrip- 
tion of the Brownian motion. The resulting equation of motion for the density operator is 
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a stochastic Liouville equation and the statistical average of the solution gives the reduced 
density matrix, the key quantity defining the system. Like the classical counterpart, the 
Langevin equation, the stochastic Liouville equation offers a convenient way for the numer- 
ical simulation of quantum Brownian dynamics, however, its efficiency is seriously limited 
due to the slow convergence of stochastic averaging 19|, |2l|, |22j. It is therefore desirable to 
derive the equation of motion for the reduced density operator or the master equation if it 
exists, given the stochastic Liouville equation. A general procedure was suggested in [l8| 
and the detailed derivation of the master equation for the dissipative harmonic oscillator 



23]. This paper tackles the dissipative dynamics of quantum optics in the 



was presented in 
same light. 

We first worked out the bath-induced stochastic fields comprising two terms with the 
rotating-wave-approximation type interaction and then showed how to determine the "dis- 
sipation operators" for a single cavity mode. Similar to the case of the dissipative harmonic 
oscillator described by the Caldeira-Leggett model, the coefficients of the master equation 
for single cavity mode are determined by a set of integral equations. It is shown that our 
result is identical to that derived by virtue of path integral technique 24j. The master 
equation of a dissipative cavity mode at a driving field was also derived and the display 
between the dissipation and the driving field was pointed out. To show that the stochastic 
formulation is a systematic method for treating dissipative dynamics in quantum optics, we 
finally explained how to acquire the master equation for the spontaneous decay of two-state 
atoms coupled to the vacuum field. For solving the master equation, because the integral 
equation is time-nonlocal, it would be better to transform it into a differential one for a 
favorable numerical implementation, if such a transformation is available. 

There are still many interesting questions in the stochastic formulation of dissipation. 
A related one to the derivation of the master equation is for what kinds of system and 
couplings such an equation exists. Notwithstanding, as the quantum dissipation becomes 
an important and subtle issue and attracts more and more attention in the community of 
quantum optics and quantum information, it is expected that the stochastic description will 
be a powerful tool in either theoretical analysis or numerical simulations. 
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